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Abstract. We investigate the dynamic critical exponent of the two-dimensional Ising 
model that is defined on a curved surface with constant negative curvature. By using 
the short-time relaxation method, we find the quantitative alteration of the dynamic 
exponent from the known value for the planar Ising model. This phenomenon is 
attributed to the fact that the Ising lattices embedded on negatively curved surfaces 
act as those in an infinite dimension, thus yielding the dynamic exponent deduced from 
the mean-field theory. We further demonstrate that the static critical exponent for the 
correlation length exhibits the mean-field exponent, which agrees with the existing 
results obtained by canonical Monte Carlo simulations. 
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1. Introduction 

Physics on curved surfaces is an intriguing subject especially with regard to second- 
order phase transition. This is because the non-zero curvature of the underlying surface 
converts the essential symmetries of the embedded physical system through a spatial 
fluctuation of the metric. The two-dimensional (2d) Ising model defined on a curved 
surface is an exemplary system. While the interest in this model originates from 
its relevance to the quantum gravity theory [U El El III E], it is presently motivated 
by the remarkable advances in nanotechnology that enables the production of curved 
magnetic layers of desired shapes [HI Ej. Earlier studies have revealed that the Ising 
model defined on a curved surface exhibits peculiar critical behaviors that differ from 
those of the planar Ising model [HI El EEH HH EH ESI- In particular, on a surface with 
a constant negative curvature, several static critical exponents have been proven to 
deviate quantitatively from the exact solution for the planar model [Hj. The latter 
results indicate the occurrence of a novel universality class of the 2d Ising model induced 
by surface curvature. 

Successful findings on the static critical properties motivate the study of the 
dynamic critical behaviors of the Ising model on a curved surface. The dynamic scaling 
hypothesis states that, in the vicinity of the critical temperature T c , the thermodynamic 
quantities O of finite-sized systems obey the dynamic scaling form as follows |15[ llfij: 



where t is the dynamic time variable, e = \T — T c \/T c the reduced temperature, L 
the linear dimension of the system, and O the universal scaling function. <f> and v are 
static critical exponents, while z is the dynamic critical exponent that takes the value 
of z ~ 2.2 in the planar Ising model (See Refs. fUEHEHj and references therein). The 
main concern is whether the finite surface curvature leads to a change in the value of z 
given above. It should be noted that the curvature-induced alteration of static critical 
exponents observed in Ref. ^1] has no bearing on this problem, since the value of z that 
characterizes the dynamic universality class of the system is generally independent of 
its static universality class. Determination of z is, therefore, crucial to obtain a better 
understanding of the curvature effect on the critical behavior of the 2c? Ising model. 

In the present study, we numerically investigate the dynamic critical exponent of 
the Ising model defined on a curved surface with constant negative curvature. The short- 
time relaxation (STR) method ^Hl EDI El E2] and the finite-size scaling analysis reveal 
that the dynamic exponent on the curved surface exhibits a value that is different from 
that for the planar Ising model. This quantitative change in the dynamic exponent 
results from the peculiar intrinsic geometry of the underlying surface on which Ising 
lattices act similar to those in an infinite dimension; therefore, the dynamic exponent 
on negatively curved surfaces yields the mean-field behavior in the thermodynamic limit. 
In addition, we observe that the static critical exponent and the critical temperature 
determined by the STR method also exhibit the mean-field behavior, as has been 
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demonstrated by canonical Monte Carlo (MC) simulations. The latter finding supports 
the conclusion of our previous work ^3] . 

2. Ising lattices on a curved surface 

In order to extract the curvature effect on the critical behavior of the Ising lattice model, 
it is desirable to adopt a simply connected surface with a constant curvature. Although 
a spherical surface is an optimal geometry, it precludes from taking the thermodynamic 
limit while maintaining its positive curvature, since it reduces to a flat plane at this 
limit. Instead, we focus on a curved surface with negative constant curvature - - a 
pseudosphere (23123- The pseudosphere is a simply connected infinite surface in which 
the Gaussian curvature at arbitrary points possesses a constant negative value. Hence, 
it is a suitable geometric surface for the consideration of the curvature effect on the 
critical properties of a system. Although the Ising models embedded on a pseudosphere 
have been considered thus far [23 123 123 I2H] , their critical dynamics have yet to be 
explored. The pseudosphere has also been considered with respect to various physical 
issues, where its intrinsic geometry is relevant to the nature of the system. These issues 
range from quantum physics |23 EDI , the string theory [SI] to cosmology jH2J ■ 

It is intriguing that the constant negative curvature of a pseudosphere permits the 
establishment of a wide variety of regular lattices [21]. The family of possible lattice 
structures consists of regular p-sided polygons that satisfy the following inequality: 



where q is the number of polygons that meet at each vertex. The series of integer sets 
{p, q] that satisfy (J2J) results in an infinite number of possible regular tessellations of a 
pseudosphere. This is in contrast to the cases of a flat plane and a spherical surface, 
where only a few kinds of regular tessellations are possible to realize. 

Amongst the infinite number of choices, we focus on a heptagonal {7, 3} tessellation 
in order to construct the Ising lattice on a pseudosphere. The resulting lattice is depicted 
in Fig. 1 in terms of the Poincare disk model, which is a compact representation of a 
pseudosphere. We emphasize that all the heptagons presented in Fig. 1 are congruent 
with respect to the metric on the disk that is expressed as follows: 



in the polar coordinate system. The metric © ensures that the Gaussian curvature 
is n = — 1 at arbitrary points on the disk (See Appendix). In addition, it leads to the 
conclusion that the circumference of the unit circle displayed in Fig. ^ corresponds to an 
infinite distance from the origin r = 0. Therefore, the regular heptagonal lattice formed 
in the Poincare disk can assume an arbitrary large system size, thereby facilitating the 
measurement of physical quantities at the thermodynamic limit while keeping curvature 
constant. 



(p-2)(g-2) >4, 



(2) 
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Figure 1. Regular heptagonal lattice established on the Poincare disk. The number 
of concentric layers of heptagons is L = 3 in this figure. All heptagons depicted within 
the circle are congruent with respect to the metric given in Eq. @ . The circumference 
corresponds to an infinite distance from the center of the circle. 



The size of our lattice is characterized by the number of concentric layers of 
heptagons, denoted by L, which effectively functions as the linear dimension in our 
lattice. It is noteworthy that the total number of sites N for our heptagonal lattices 
exhibit a complicated dependence on L that is expressed as 

N(l) = 7, 

L-2 

N(L) = 7 + 7j2 (u+v+ j + U-vJ) , [L > 2] (4) 
where u± and v± are defined by 

u± = 2± v / 5, v ± = ^±. (5) 

When L ^> 1, the total sites N grows asymptotically as N oc e L , which is quite rapid 
in comparison with the case of the planar Ising model: iV oc L 2 . This exponential 
increase in N is a manifestation of the constant negative curvature of the underlying 
geometry, resulting in that the ratio [N(L) — N(L — 1)]/N(L) approaches a non-zero 
constant 1 — e _1 in the limit L — > oo. The latter feature means that the boundary sites 
of our lattices can not be neglected even in the thermodynamic limit, but contain a 
finite fraction of the total sites. 

Boundary effects coming from these sites are difficult to be eliminated completely, 
because the periodic boundary conditions are hard to be employed to regular lattices 
assigned on a pseudosphere. Thereby, we have used the following manner in order to 
reduce the boundary effects on physical quantities of the system. Suppose that an Ising 
lattice consists of L ont concentric layers of heptagons. Then, for employing the STR 
method (see Sec. 3 for its details), we take into account only the Ising spins involved in 
the interior L in layers (L in < L out ) so as to reduce the contribution of the spins locating 
near the boundary. In actual calculations, L m is varied from 4 to 6, and for each L m the 
number of disregarded layers 5L = L out — L m is systematically increased from to 3. 
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By examining the asymptotic behavior for large 5L, we can deduce the bulk properties 
of the Ising lattice model embedded on the pseudosphere. Similar procedure has been 
employed thus far [213 ES| to extract the bulk critical properties of the negatively curved 
system. 

3. Short-time relaxation method 

Our aim is to investigate the dynamic critical behavior of the ferromagnetic Ising model 
defined on the heptagonal lattice. The Hamiltonian of the system is given by 

H — — J ^2 s * s i> s i = ±:L > ( 6 ) 

<i,j> 

where denotes a pair of nearest-neighbour sites and J(> 0) is the coupling strength 
between them. Hereafter, temperatures and energies are expressed in units of J / k-Q and 
J, respectively. 

In order to determine the dynamic critical exponent, we use the short-time 
relaxation (STR) method [TBI E01 EH l22j • This method enables neglecting equilibration, 
thereby avoiding computational burden significantly as compared to dynamical 
simulations in equilibrium. The basis of the STR method is the measurement of the 
quantity [THJ ESI 

Q(t) = (sgn ^ , sgn(x) = | _ j j* > °j« (7) 

which is obtained by the conventional MC procedure. Here, the time t is measured in 
terms of one MC sweep, and the angular bracket indicates that the average is taken 
over different time sequences that begin from the same initial configuration. In actual 
simulations, the initial configuration is selected as Si = 1 for all sites, thereby yielding 
Q(0) = 1. It then follows a monotonic decay with an increase in t, and finally Q(oo) = 0, 
since there is no preferred direction at equilibrium (t — ► oo). 

A dynamic scaling hypothesis suggests that in the vicinity of the critical 
temperature T c , the quantity Q obeys the following scaling form: 

Q(t, e,N) = Q (tN~ g , eiV 1 ^) . (8) 

Here, z and /i are referred to as the dynamic and static critical exponent, respectively, 
for the Ising lattices on the pseudosphere (The physical meaning of these exponents will 
be given soon below). It is noted that N instead of L is employed in (JBJ) as a scaling 
variable; this is a direct consequence of the exponential increase N oc e L discussed in 
Sec. 2%. 

The dynamic exponent z is determined as follows. At the transition point T = T c , 
e = so that the second argument eN 1 ^ of the scaling function Q vanishes. This yields 

Q(t,N)\ e=0 = f(x) (9) 

| Similar argument has been made regarding an infinitely coordinated Ising model j.'i.'ij and a randomly- 
connected Ising model [21, where N is the only variable determining the system size. 
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with the definition x = tN~ z . Equation (jUJ) states that, if we calculate Q(t) by 
maintaining T = T c and plot the results against x by selecting an appropriate value 
of z, the curves of Q{t) for different iVs should collapse onto a single curve. Hence, the 
dynamical exponent z can be identified as the optimal fitting parameter for Eq. ©. We 
should note that the evaluation of z by the above procedure requires the knowledge of 
the value of T c beforehand; therefore, in actual calculations we use the numerical data 
of T c obtained from canonical Monte Carlo simulations for the same model [T4*] . 

Apart from the value of z, the STR method enables us to determine the static 
critical exponent //. This is achieved by extracting the data Q(to) at the specific time 
t = to for various T and N. The time to is defined such that parameter a = toN~ z with 
fixed z is invariant to the change in N. Under this condition, we obtain 

Q^N)\ t=tQ =g(y) : (10) 

where y is defined as y = eN 1 ^. As a result, fitting the data of Q(to) for different T 
and N onto a single curve against y yields the values of /i and T c as the optimal fitting 
parameters. This procedure provides a complementary method to determine \x and T c , 
which were extracted in [T4*] . 

It deserves further comment about the exponents z and \x introduced in (jHJ). From 
the scaling hypothesis, the exponent z is assumed to describe the dynamic scaling law 

rocfr* (11) 

between the relaxation time r and the correlation volume £y. Here, £y is a natural 
generalization [Hill EE] of the localization length £, thus obeying the power law £y oc e _At 
in the vicinity of T c . When the underlying geometry is a <i-dimensional flat surface (or 
space), £y oc £ d holds so that (JTTJ) yields z = z/d ~ 1.1 for two-dimensional planar 
Ising models§. However, when the underlying geometry is curved, the relation £y oc £ d 
becomes invalid so that z may take a value different from the above. Similarly, the static 
exponent \i reduces to — vd — 2 for planar Ising models (because v — 1 |HE]), while it 
is not the case for curved cases. In fact, our numerical simulations have revealed that 
the dynamic exponent z for heptagonal Ising lattices exhibits a quantitative deviation 
from that for planar ones, as will be seen in the next section. 

4. Results: Dynamic exponent z 

4-1. z for entire lattices (SL = 0) 

Before addressing the bulk critical properties, we first demonstrate the results for entire 
heptagonal lattices with 5L = (boundary contributions are fully involved). Figure |2] 
(a) shows the short-time relaxation behavior of Q(t) for various systems sizes N. The 
temperature is fixed at T = T c = 1.25 consistent with Ref. and the standard local 
update algorithm |HZ| is employed in order to calculate the time evolution of Q(t). N 

§ In d-dimensional planar Ising models, £y oc £ d so that r oc £ z = £y z / d = £y z , and £y oc = e~ vd 
since £ oc (r v . Eventually we obtain z — z/d and y, = vd. 
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t [10 3 MCS] x = tN ' z Disregarded layers: AL 

Figure 2. (a) Calculated results of Q(t) for various lattice sizes N. (b) Finite size 
scaling of Q(t,N) against the scaling variable x = tN~ z . The critical temperature 
T c = 1.25 is assumed in accordance with Ref. |14|. The estimated value of z is 
z = 0.90(2) with a 95% confidence interval, (c) The dependence of z on the number 
of disregarded layers SL. The solid line serves as a guide to eye. 

is varied from 315 to 2240; this corresponds to the change in the number of concentric 
layers from L in = 4 to L m = 6 (See Eq. (HJ)). For each value of N, we have performed a 
sample average for more than 3 x 10 4 independent runs. 

We see that all curves of Q(t) in Fig. El (a) exhibit a monotonic decay from the 
initial value of Q(0) = 1 toward the equilibrium value Q(oo) = 0. The decay time grows 
systematically with an increase in the system size N as expected from the scaling law 
r oc N z . Hence, rescaling the horizontal axis t of Fig. El (a) by dividing it by the factor 
N z with an appropriate exponent z yields the fitting of all curves into a single curve. 

Figure El (b) presents the scaling attempt of Q by taking x = tN~ z as a scaling 
argument. The dynamic exponent z is determined as the optimal value that minimizes 
the following quantity [TTHj : 

(12) 

where a, and \ are fitting parameters, and o~(x) is the standard deviation of Q. Note 
that in (|12|). the entire x range is partitioned into n seg segments Xi with equal intervals. 
We have set n seg = 20 in actual calculations, and extracted the values of and i>; 
(1 < i < 20) that minimize A(z) for a given z. By probing the behavior of A(z) with 
respect to z, we have attained the optimal value z = 0.90(2) with a 95% confidence 
interval. The resulting plot shown in Fig. El (b) displays the best collapse onto a single 
curve f(x) in a broad range of the scaling variable x = tN~ z . Importantly, the resulting 
value z ~ 0.9 is somewhat smaller than that of the planar Ising model z — z/d ~ 1.1. 
Since the boundary contribution is fully involved in the present case (SL = 0), the above 
deviation of z possibly stems from the mixture effect of the boundary Ising spins and 
the interior heptagonal Ising lattice. 



Mz) 



EE 

i=l XGX; 



Q(x) - di - hx 



o~[x ) 



The dynamic exponent of the Ising model on negatively curved surfaces 



8 



4-2. z for bulk lattices (SL > 1) 

We now turn to the study of bulk critical properties of the heptagonal Ising lattice 
model. As mentioned in Sec. 2, the boundary spins of Ising lattices on a pseudosphere 
are thought to affect significantly to the nature of the system; this is because the 
number of the spins along the boundary increases as fast as that of total spins of the 
lattice. Thereby, in order to extract the bulk critical exponents, we try to remove the 
boundary contribution by setting the disregarded layers SL = L ont — L UI to be finite, 
i.e., by summing up only the spins within the interior L in layers when performing MC 
simulations on the systems consisting of L out (> L in ) layers. 

The calculated results are summarized in Figure 0(c); each data of z was extracted 
by the scaling analysis for lattices of 4 < L in < 6 and a given SL. The appropriate 
value of T c for each SL is again referred to the results obtained from canonical MC 
simulations [14J. (Hence, the maximum system size we have treated reaches L ont = 9, 
which corresponds to N = 40432.) We have observed that z monotonically decreases 
with increasing SL within the range < SL < 3. The monotonic decrease in z indicates 
that for sufficiently large SL, z takes a value that is considerably smaller than the 
corresponding value of the planar Ising models: z ~ 1.1. Thus, it follows that the bulk 
Ising lattices assigned on the negatively curved surfaces belong to a dynamic universality 
class distinct from those of the planar Ising models. The determination of the asymptotic 
value of z for 5L ^> 1 would provide further conclusive information; this would require a 
huge computational effort. Instead, we have provided a concise argument regarding the 
asymptotic value of z in Section 6; this argument implies that the dynamic exponent in 
our system takes the mean-field exponent ^mf = 1/2 in the limit SL — > oo. 

5. Results: Static critical exponent // 

We now focus on the estimation of the static critical exponent /i and the critical 
temperature T c on the basis of the scaling relation (jlOjl . Figure El (a) presents the 
T-dependence of Q(T) under the condition SL = 0. Parameter a = tN~ z is fixed at 
a = 1.5, where z = 0.90 is set in accordance with the result from Fig. 2 (b). We observe 
a unique crossing point at Q ~ 0.74 and T ~ 1.25. Even when the values of parameter a 
are varied, a crossing point still appears at a temperature almost identical to T ~ 1.25, 
as expected from Eq. (JSJ). Hence, it provides a rough estimate of the critical temperature 



The full scaling plot for Q(T) is displayed in Fig. El (b). The fitting process was 
based on a polynomial expansion of the scaling function g(y) expressed in Eq. (fTU|) as 



where Cj, T c , and fi are considered as fitting parameters. A very smooth collapse is 
obtained with /x = 3.45(4) for the upper branch, and \x = 3.36(3) for the lower one. 
The critical temperature is estimated as T c = 1.254(3) for both branches. It was also 
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Figure 3. (a) Calculated results of Q{T) for various values of N. Parameter a — tN~ z 
and the exponent z are fixed at a = 1.5 and z = 0.90 for all values of N, respectively, 
(b) A scaling plot of Q(T,N) against the scaling variable y = eN 1 ^. The optimal 
values of fx for upper and lower branch are Li = 3.45(4) and fx = 3.36(3), respectively. 
The critical temperature is estimated as T c = 1.254(3) for the two branches, (c) 
<5L-dependences of /i and T c . 



confirmed that the values of \i and T c are almost invariant to the change in the value of 
parameter a that is arbitrarily chosen within the range of 0.5 < a < 3.0. 

We have then clarified the boundary contribution to the values of fi and T c by 
setting 8L > 1. Figure El (c) shows the <5L-dependence of v and T c under the same 
numerical conditions with respect to L in and a as in the case of 5L = 0. It is found 
that both /i and T c decreases monotonically with increasing 5L. Further noteworthy is 
the fact that all the values of fi and T c given in Fig. El (c) are in excellent agreement 
with those obtained by MC simulations in equilibrium [J jJ. where the tendency of [i 
toward the mean field exponent /imf = 2 in the limit 5L — > oo has been suggested. This 
consistency with respect to /i and T c strongly supports the validity of our STR approach 
to the critical behavior of the Ising model on curved surfaces. 

6. Discussions 

It has been demonstrated that the dynamic critical exponent z of the heptagonal Ising 
lattice model exhibit a quantitative deviation from that of the planar Ising model 
z ~ 1.1. In fact, it takes a value ~ 0.9 under the condition 5L = 0, and still continues to 
decrease almost linearly with 5L. This behavior of z indicates the occurrence of a novel 
universality class with regard to dynamic criticality, and evidences the conjecture that 
the peculiar intrinsic geometry of the underlying surface affects the dynamic critical 
properties of the embedded system. Besides, the static critical exponent n also turns 
out to exhibit a monotonic decrease with increasing 5L toward a particular value of 
Hmf — 2, which is consistent with the value derived from the mean field theory. 

The mean field behavior of fi has been suggested earlier by Rietman et al. j2H] 
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and by Doyon and Fonseca |2Ej- This behavior is attributed to the fact that an Ising 
lattice embedded on a pseudosphere is effectively an infinite dimensional lattice at large 
distance due to the exponential growth of the total spins iV [39]. For an ordinary Ising 
lattice in d dimension, the number of spins along the boundary, N B , is expressed by 
N s oc jV 1_ ( 1//d ). Hence, by comparing it with the peculiar relation N s (x N that is 
valid on negatively curved surfaces, we eventually reach the consequence d = oo for the 
latter geometry. As a result, our heptagonal Ising lattices yield the mean-field critical 
exponent /zmf = = 2 [SSI, where z/mf = 1/2 is the mean-field exponent and d c = 4 

is the upper critical dimension for the Ising model. 

In this context, the dynamic critical exponent z on negatively curved systems is 
supposed to exhibit the mean field exponent z MF = z M f/4 = 1/2 at 5L ^> 1, since 
zmf = 2 for the Ising model. This inference is indeed consistent with our numerical 
results; as shown in Fig. 2 (c), z takes a smaller value than the planar case (z ~ 1.1) 
whatever value 5L is, and moreover, it monotonically decreases with reducing the 
boundary contribution up to 5L = 3. Hence, when extrapolated to the bulk limit 
5L ^> 1, it yields the mean field exponent zmf — 1/2 as similar to the case of static 
critical exponents [T3j. We also note that it is interesting to investigate the dynamics 
exponents associated with the Griffiths-type phase transition in the Ising model defined 
on negatively curved surfaces [27] ; a comparison between the dynamics exponent of the 
bulk region and that of the boundary region will provide a better understanding of the 
critical dynamics of entire Ising lattices defined on negatively curved surfaces. 

Further noteworthy is that the surface curvature effect on critical properties of the 
assigned spin lattice model would become richer if the spin variables possess orientational 
degrees of freedom. In the latter cases, the interacting energy of neighboring spins is 
determined by their relative angle, which is a function of a spatially-dependent metric 
tensor. Thereby, the Hamiltonian of the system involves the metric tensor in an explicit 
form. As a result, the energetically preferable spin configurations and the temporal 
fluctuations of the order parameters are directly influenced by the intrinsic geometry of 
the surface determined by the metric tensor. It is thus expected that these vector-spin 
models defined on curved surfaces exhibit peculiarities associated with phase transition 
and low energy excitations (101 EH E21 EHl H3]. To the best of our knowledge, there 
was no attempt to reveal the dynamic critical exponent of vector-spin lattice models 
defined on curved surfaces. The STR method that we have utilized is applicable to these 
systems; detailed analyses of this issue will be presented in future. 

7. Conclusion 

In the present study, we have investigated the dynamic critical behavior of the regular 
heptagonal Ising model defined on a curved surface with constant negative curvature. 
The STR methods that incorporate finite-size scaling analyses have been employed in 
order to compute the dynamical critical exponent z as well as the static critical exponent 
/j, for the correlation volume £y. The resulting values of both the exponents are distinct 
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from those of the planar Ising model. Furthermore, we have studied quantitatively how 
the boundary spins contribute to the determination of z and //, and eventually we have 
revealed that they both reduce to the mean field exponent zmf = 1/2 and /xmf = 2 in 
the bulk limit. We hope that our results prove to be a fundamental basis for further 
fruitful studies on the critical behaviors occurring on curved surfaces. 
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Appendix A. The Gaussian curvature in the Poincare disk 

This appendix demonstrates that the Poincare disk possesses the negative constant value 
k — — 1 at an arbitrary point within the disk. The Gaussian curvature k at a given 
position is defined by jlHj 

K = -R/2, (A.l) 

with the quantities: 

R = g lj Rij, Rij = g kl R ki ij . (A. 2) 

Here, R, R^, and Rkiij are termed as the scalar curvature, Ricci tensor, and curvature 
tensor, respectively. (Hereafter, we will use the summation convention on repeated 
indices.) The g^ is the inverse of the metric tensor that determines the line element 
at a given position as ds 2 = gijdx % dxK Note that all the three quantities R, R^, and 
Rkiij are functions of gij. Hence, in principle, the knowledge of and the explicit form 
of Rkiij are sufficient to calculate the value of k, while the actual calculations are lengthy 
as shown below. 

The explicit dependence of the curvature tensor R k uj on g^ is expressed as: 

Rkiij = 9kmK"ilji (A. 3) 

R m Uj = d l r™ i - d.n + r m pl r% - r™ pj r% (aa) 

where i" 1 ™- is the connection defined by 

mp 

r % = — (d j9pi + d t g PJ - d p9ij ) . (A.5) 



Equations ()A.3|) - ()A.5|) yield an alternative form of Rkiij, which is expressed as: 

R m = l(dA9M + d k d igij -d k d jgu -dA 9kj ) 

+ 9rs (r lk r% - r Jk r u ) . (A.e) 

It is noteworthy that due to the symmetric property of Eq. (jA.6|) . the tensor R k uj has 
only four non-zero components in a two dimensional system. Further, these components 
are related to each other in the following manner: 

-^1212 = —RlXVl — —R\22\ — -^2121- (A-'' 7 ) 

As a result, the tensor Rkuj is determined by only a single component Rui2- in addition, 
the symmetric property of Rkiij results in the following equality with respect to the 
components of the Ricci tensor R^: 

Rn = g 22 R2i2i, R22 = g n Rm2, (A. 8) 

and Ru = R21 = 0. Consequently, we can represent the scalar curvature R in two- 
dimensional systems as 

R = g n R n + g 22 R 22 



g g -K2121 + g g ^1212 

2-Rl212 



(A.9) 



The dynamic exponent of the Ising model on negatively curved surfaces 



13 



where g = det^], and thus g 11 = (722/(7 an d g 22 = gn/g- Equation (jA.9|) yields the key 
relation between the component -R2121 and the Gaussian curvature k as follows: 

k = -^. (A.10) 
9 

This identity is known to hold for general two-dimensional systems 

The value of -R1212 for our model is evaluated by the following procedure. For the 
Poincare disk, the metric tensor gij is represented in the matrix form as follows: 

/ 4/ \ 1 
M = . = 77—^- (A.11) 



4r 2 / 



1 - r 



This means that the line element on the Poincare disk is given by ds 2 = 4f(dr 2 + r 2 d8 2 ), 
as already shown in Eq. (jHJ). Straightforward calculation based on the metric (jA.llJ) 
reveals that there are four non-zero components of the Christoffel symbol expressed 
as 

g u g 11 
r\i = -y^i^i, r\ 2 = — —8^22, 

22 

r\ 2 = r 2 21 = 9 —d l9 22. (A.12) 

By substituting Eqs. (jA.llJ) and ()A. 12j) into the expression of the tensor Rkiij (|A.6|) . we 
directly obtain: 

16r 2 1 

-R1212 = -r, 2T4 = ~- (A.13) 

(1-r 2 ) 4 g 

Eventually, by comparing it with Eq. (jA.lOJ) . we arrive at the final conclusion k — — 1 
for arbitrary positions within the Poincare disk. 
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